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1.  Introduction 

This  paper  presents  a  mathematical  theory  for  the  deforma- 
tions of  isotropic  shells  under  the  applications  of  edge  forces. 
The  theory  is  obtained  by  imposing  the  principal  of  virtual  work 
on  the  class  of  deformations  which  satisfy  the  Kirchhoff 
hypothesis,  i.e.  those  deformations  which  carry  normals  to  the 
undeformed  middle  surface  into  normals  to  the  deformed  middle 
surface  with  no  change  in  length  along  normals.   It  is  not 
assumed  that  displacements,  strains,  or  slopes  are  small. 
Although  in  the  derivation  of  the  model  it  is  not  assumed  that 
the  shell  is  thin,  v;e  do  not  expect  the  theory  to  be  physically 
realistic  unless  the  shell  is  in  fact  thin. 

Imposing  the  principal  of  virtual  work  yields  a  system  of 
differential  equations  for  the  three  components  of  position  of 
the  deformed  middle  surface  and  six  relations  between  the 
applied  edge  forces  and  the  deformed  middle  surface.   It  will 
be  shox-m  that  the  system  of  differential  equations  can  be 
expressed  as  a  tenth  order  system.   The  six  relations  between 
edge  forces  and  deformed  middle  surface  contain  two  arbitrary 
functions  and  hence  represent  four  boundary  constraints  on  the 
deformed  middle  surface  when  the  edge  forces  are  prescribed. 
It  will  be  shown  that  to  be  physically  realistic  one  of  the 
arbitrary  functions  should  be  chosen  to  be  zero.   We  then  have 
five  boundary  constraints  for  the  tenth  order  system  of 
differential  equations.   Hence  after  prescribing  the  edge  forces 
we  expect  the  deformation  to  be  uniquely  determined  unless 
buckling  occurs. 
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Since  it  has  recently  been  proved  rigorously  [l]  that  the 
nonlinear  equations  of  the  three  dimensional  theory  reduce  under 
suitable  hypotheses  to  the  equations  of  von  Karman  and  P8ppl 
for  plates,  it  seems  essential  that  any  nonlinear  shell  theory 
which  purports  to  be  more  general  than  the  von  Karman  and  F8ppl 
theory  should  also  reduce  to  that  theory  when  one  makes  suitable 
restrictions.   Unfortunately,  as  will  be  shown  in  this  paper, 
the  theory  derived  on  the  basis  of  the  hypothesis  stated  above, 
does  not  meet  this  test.   It  is  felt  that  it  is  worthwhile  to 
point  this  out.   In  fact,  Naghdi  and  Nordgren  [2],  for  example, 
have  recently  presented  a  theory  for  the  treatment  of  thin 
shells  in  which  the  essential  hypothesis  is  the  same  as  that 
of  this  paper.   It  is  not  easy  to  decide  directly  whether  the 
the  theory  of  Naghdi  and  Nordgren  is  equivalent  to  the  one 
presented  in  this  paper,  but  the  author  has  checked  the  point 
regarded  here  as  important,  i.e.  that  the  theory  of  Naghdi  and 
Nordgren  does  not  yield  the  von  Karman-P5ppl  equations  with  the 
correct  coefficients  when  it  is  formulated  for  thin  plates. 
Thus  it  also  does  not  agree  with  the  von  Karman-F8ppl  equations. 
The  fact  that  others  are  continuing  work  on  theories  of  this 
type  increases  the  importance  of  calling  attention  to  this 
defect  in  such  theories. 

2.  The  Three  Dimensional  Theory 

In  deriving  the  shell  theory  we  attempt  to  copy  a  deriva- 
tion of  the  complete  three  dimensional  nonlinear  theory.   In 
this  section  vje  present  that  derivation  of  the  three  dimensional 
theory. 


n^. 


Consider  a  fixed  rectangular  Cartesian  reference  frame  X. 
Let  (x, jXojX^)  be  the  position  of  a  particle  in  the  undeformed 
state  and  (x^^XpjX.,)  be  the  position  of  the  same  particle  in  the 
deformed  state.   The  Eulerian  coordinates  x.  are  considered  as 

functions  of  the  Lagrangian  coordinates  x.. 

Sx. 
Let  p.  .  =  -rr^T^  and  let  W  =  W(p.  .,x.  )  be  the  strain  energy 

per  unit  undeformed  volume.   The  elastic  properties  of  the 

material  are  then  determined  by  the  choice  of  W. 


Let  E  =   j   VJdv  be  the  total  strain  energy  where  R  is  the 

r' 

undeformed  region  and  dv  is  the  element  of  undeformed  volume. 

Let  (t-,,tp,t^)  be  the  force  vector  applied  to  the  boundary 
S  of  the  deformed  region  R.   The  units  of  t.  are  force  per  unit 

deformed  area. 

rr     -    - 
The  principal  of  virtual  work  states  that  E  =  //  t.x.dS 

S 

where  the  summation  convention  is  used,  a  super  dot  denotes  the 

variation,  and  dS  is  the  element  of  deformed  area. 

c^i  •     rrr       ^^i 

Letting  q.  -  =  -r-^ — ,  we  have  E  =  ///  q.  .  v —  dv.   Using  the 
ij   '^y.:  ,•  JJ'^J      ij  ex,- 


ij  V        J 


divergence  theorem  this  becomes  E  =  //  x.q.  .n.dS-  ///  x.  ^  ■'*''-'  dv 

J J       1^1 J  J    JJJ       1  ox. 

S  R        "J 

where  S  is  the  boundary  of  R,  (n^,np,n^)  is  the  unit  outer 

normal  vector  to  S,  and  dS  is  the  element  of  undeformed  area. 

Hence  the  principal  of  virtual  work  yields  the  differential 

equations 

hq.  . 

(2.1)  -..-iJ-  =  0 

ox  . 

3 

for  the  functions  x.  and  the  relations 

1 
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(2.2)  t.  =  q.  .n.  -^ 


between  the  force  and  deformation. 


Letting  {n^,n2,n,)  be  the  unit  outer  normal  vector  to  S, 
it  can  be  shovm  that  q.  .n.  —  =  -. — p  q.  p,  .n  vjhere  |p|  is  the 
determinant  of  the  matrix  p  =  (p.,).   Thus  if  we  let 

'^•''  -^ij  =  TFT  ""'^j!^ 

the  relations  between  force  and  deformation  become 

(2.4)  t.  =  T.  ,n.  . 
Also  one  can  show  that 

OT  .  . 

(2.5)  ^ii  =  0  , 

and 

(2.6)  T.  ,  =  T  ,, 

the  latter  being  true  if  VJ  is  invariant  under  rigid  displacements. 
Thus  one  obtains  the  usual  differential  equations  and  symmetry 
for  the  components,  t.,.,  of  the  stress  tensor  and  also  the  usual 
relationship  bet\jeen  the  t.  .'s  and  the  applied  forces. 
For  a  homogeneous  isotropic  material  we  have  [5] 

(2.7)  q^j  =  Vio  +  Vio^  Vik^^kP^J  • 

Here  the  W. 's  are  functions  of  the  principal  strains  and  the 

c^  .'s  are  the  components  of  the  rotation  matrix  c.   Specifically 

by  principal  strains  vje  mean  the  three  stationary  values  of 


en: 


•;ii5 


ds  — 

^  -  1  where  ds  and  ds  are  arc  length  in  the  undeformed  and 

deformed  bodies  respectively.   Letting  s.  be  the  sum  of  the  1-th 

poivers  of  the  principal  strains,  v;e  have  W,  =  2  !  -^ —  -  3  ^—), 

,,    ci;    „  d\l      ,  ^  SW  ,,  ^  2'J    _     ^   ^'^      '^     ^    .      ^' 

V/„  =  ^^ —  -  2  -r —  +  J>  -^ — ,  W_  =  3  ^; —  .  The  rotation  matrix 
2   ds,      cso     ^s-2        3     OS-. 
1.2       3  ;; 

c  is  defined  by  requiring  that  c  be  orthogonal,  det  c  =  1,  and 
c*p  (or  equivalently  pc*)  is  symmetric  and  positive  definite. 
The  matrix  c  then  rotates  the  directions  of  principal  strain 
in  the  undeformed  body  into  the  corresponding  directions  of 
principal  strain  in  the  deformed  body. 


Also  c*p  =  j p *p  so  that  the  eigenvectors  of  c*p  have  the 

direction  of  principal  strain  in  the  undeformed  body  and  its 

d's 
eigenvalues  are  the  stationary  values  of  -r-  . 

Cl  s 

In  terms  of  p  and  c,  s,  =  [(c*p-l)"'"]  where  the  square 

1 

bracket  denotes  the  trace  of  the  matrix. 

The  above  expression  for  q.  ,  will  be  used  in  deriving  our 
theory  for  isotropic  shells. 


3»  Terminology  for  Shell  Theory 

Let  r  =  v{9     9    )   be  the  position  vector  to  the  undeformed 
mj-ddle  surface,  9     and  9^   being  parameters  for  the  middle 
surface.   Let 

(>1)       Sa  =  ^,c.   ^"^   f 3  =  --— zr-  • 

Igl^  ggl 

Here  and  in  all  that  follows  Greek  indices  take  on  the  values 
1,  2  and  Latin  indices  take  on  the  values  1,  2,  3.   Also 
,a  =  -rr^   and  ,i  =  -'-g— •  Observe  that  g_  is  a  unit  normal  vector 
to  the  undeformed  middle  surface. 


:>  Of;  T 
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We  Introduce  curvilinear  coordinates  9    ,    9        9      in  the 

±  d  ^ 

undeformed  shell  by  letting 

(3.2)  3?=  (x^.XgjX^)  =  r+9^     , 

Obcerve  that  r  and  g'.  are  all  independent  of  0^.      Then  9^ 
is  the  distance  from  (x-,^Xp,x-,)  to  the  undeformed  middle  surface 
along  a  normal  to  the  undeformed  middle  surface.   We  assume  that 
the  shell  is  thin  enough  so  that  normals  to  the  undeformed 
middle  surface  do  not  intersect  in  the  shell. 

Let 

(3.3)  g^j  "  ^i'^'j  '    ^s'"^)  =  ^Sij^'"^  '        ^"^"^       "S ^  "  S^'^'^j  • 

Then  the  go's  and  g   ' s  are  the  covariant  and  contravariant 

components  of  the  usual  metric  tensor  for  the  undeformed  middle 

3a  ^5 

surface.   Also  g.,  =  g   =0  and  g^^  =  g'   =  1.   Finally 

g  ••  =  3  •  g  ,  g^  •  g  J  =  5^ ,  and  g^  =  g  ^  . 

All  vectors  associated  v;ith  the  shell  v/ill  be  expressed 
in  terms  of  thej.r  components  wj  th  respect  to  the  vectors  g. 
or  g  .  This  is  done  because  we  v;lll  be  integrating  the 
components  of  vectors  with  respect  to  9...      Since  the  vectors  g. 
and  g  ""■  are  independent  of  9        such  integration  yields  the 
components  with  respect  to  the  vectors  g..  or  g  of  the  integral 
of  the  associated  vector. 

lie  express  the  vectors  g.  and  g"""  in  terms  of  our  fixed 
Cartesian  coordinate  system  as  follows: 

1  o2  .3> 


(3.*+)     t^  =   (a^,a^,ar')  ,   f ^  =  (A^,A^,A^)  . 


si'.i 


■gfic;  c:-l3?y 


ic   ■  ,J 


i59:f.r'»  "yj'..-  'ic  amic-t  nz 


.,*;-.:■A^  = 


,■  t^ 


Then 

Given  any  indexed  set  of  functions  associated  with  the 
Cartesian  coordinates  (such  as  the  p.  's  of  Section  2)  we 
associate  inde::ecl  sets  vjith  the  curvilinear  coordinates  as 
follows:   P.  .  =  a^a"p   ,  P."^'  =  a^^A'^p   ,  P^  .  =  A^a"p   , 

P-^-^*  =  A^A-^p.^  .   We  then  have  P.  .  =  g.,P^V,  etc.,  so  that  the 
m  n'-  vail  ij   '^ik  j '     •* 

quantities  g"^'^  and  g.,  .  can  be  used  to  raise  and  lower  indices 
in  the  usual  way.   Actually  this  means  the  indices  1,  2  are 
raised  and  lowered  using  the  quantities  g   and  g  ^  and  the 
index  3  has  the  same  meaning  regardless  of  whether  it  is  up  or 
down.   This  procedure  for  generating  indexed  sets  associated 
with  the  curvilinear  coordinates  and  for  raising  and  lowering 
the  indices  applies  regardless  of  the  number  of  indices  involved. 
Under  changes  of  curvilinear  coordinates  of  the  middle  surface 
v;e  notice  that  the  quantities  P  o  transform  like  the  components 
of  a  covariant  tensor  of  rank  two,  the  quantities  P„^  and  P_^ 
like  the  ccmoonents  of  a  covariant  tensor  of  rank  one,  and  P^^ 
is  invariant.   Similar  statements  hold  for  indexed  sets  with 
any  number  of  indices. 

In  Table  I  we  have  listed  in  the  left  column  certain 
indexed  sets  associated  vjith  the  fixed  Cartesian  coordinate 
system  and  in  the  right  hand  column  one  of  the  corresponding 
indexed  sets  associated  with  the  curvilinear  coordinates. 


'  ;.  .  i . 
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r;rl,J. 
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Table  I 

x-components      9-components 


Pi-,- 

P.  . 

c.  . 

=ij 

^ij 

«« 

®ap 

^aP 

®ijk 

^Ijk 

Here  e^g  =  ±1  (®i ^v  =  *1)  according  as  aP(ijk)  is  an  even 
or  odd  permutation  of  the  integers  1,  2  (1,2,5).   If  two 

indices  are  equal  e^^  =  e.  .,  =  0.   It  can  then  be  shown  that 
'■  af3    ijk 

r  -  n  ap  _  1 

^ap  "  ^  S  e^P   ,    e    -  --  e^^ 


/i 


(3.5)     ' 


.     ^ijk   ^-  ®ijk  '        ^  J-   ®ijk 

v;here  g  =  det  (g.  , ) . 

In  terms  of  this  notation  it  can  be  shown  that 


(5.6)  { 


Si^  Sij    e^jj^ 


e 


\     g  X  g  "^  =  e  '^   g 


'k 
It  is  easily  verified  that 


(5.7)        Q  '^  =  7T5 —  >   Q  i  =  T  >  etc., 

and  from  (2.7) 

(5.5)        Q.  .  =  VLP.  .  +  W„C.  .+l/^P.,C.^  P^, 
-^  ij    1  ij   2  1,1    5  ik  i    J 


for  an  isotropic  material. 


t:<; 


. .  'JO  V    ^ 


C. .  is  determined  by  the  conditions  C,  . C  .  =  g. ., 
det  (C.  .)  >  0,  and  (C,  .P^.)  (or  equivalently  (P.,C,.^))  is 

X,J  KIJ  IKJ 

symmetric  and  positive  definite. 

If  d  g.  and  D  g.  are  unit  vectors  in  corresponding 

directions  of  principal  strain  in  the  strained  and  unstrained 

bodies,  then  C-^-^D.  =  d^. 

J 

We  also  observe  the  following.  Let  £.  .  be  an  indexed  set 
associated  with  Cartesian  coordinates  such  that  H.  .  =  &...  Let 
(y,  ,yp,y-2)  be  an  eigenvector  of  the  matrix  {&.  .}  and  A  a  corre- 
sponding eigenvalue.  Let  L. .  and  Y.  be  indexed  sets  associated 
with  the  curvilinear  coordinates  derived  from  ^3.  .  and  y. 

X  J  X 

respectively.   Then  (L,.  .)  and  (L"'"'^)  are  symmetric  and  L.*^  =  L  .. 

X  ^1  X  J 

We  write  L^  for  L^  .  and  L.K      Then  L^Y-^  =  AY^  so  that  A  is  an 
eigenvalue  of  L^  and  (y-j^jyg^y^)  =  Y^g^. 

It  follows  that  if  ^'^'Ea     is  a  vector  in  a  direction  of 
principal  strain  in  the  undeformed  body  and  A  is  the  corre- 
sponding value  of  ~  ,  then  C^^P^  D-^  =  AD^  and  P^^Pj^j  D^  =  aV. 

Next  we  introduce 

the  components  of  the  second  fundamental  form  of  the  undeformed 
middle  surface,  and  vie   define  [  '..':',  by  the  equations 

It  can  be  shovm  that  the  Pp,^  s  are  the  usual  Christoffel 
symbols  for  the  metric  of  the  undeformed  middle  surface. 


j.:.n. 


[0     iBU&i: 


,0    --   .  . 
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Since  'g  -"g  •  =  5"^,  we  have 

J         J 


,a'Sj 


— ^i  p\  k  -* 

-g  •  I  Jet  gj^ 


-r ' 


and  hence 


(3.11) 


^  ,  a   ~  '  j  a  '^ 


If  u"*"  g*.  is  any  vector  depending  only  on  9,  and  9^,    we 


define  u"""  I   and  u.  |   by  the  relations 

Ct  X     Co 


Then 


/    i-*    \  i  I    — *  I    -*1 


(3.12)       u^l^  =  u\„+  TjiuJ   ,       uj^  =  u.  _„  -  rjuj , 


u    1^  =   S   ^u 


J 'a 


We  observe  that  u  1   is  not  an  ordinary  covariant 
derivative  since  the  f'^-i's  are  not  all  Christoffel  symbols 

and  the  u  's  are  not  the  components  of  a  tensor. 

ict 
For  doubly  indexed  sets  D   depending  only  on  0,  and  9p 

where  the  first  index  can  have  the  values  1,  2,    3  and  the  second 

index  the  values  1,  2,    we  find  it  convenient  to  make  the 

following  definitions. 


(3.13) 


fDi-L  =Di%-.r>^%r  «Di^ 


P 


i  v.. 


I   =  D  ""     r  -^D  ^  +  r  ^'d  ' 
'P   ^i  ,3   '  iP^j  +"  7P°i 


D^ 


a,P 


■^r.jDJ  -riD 


7T.i 


jP 


J 


a(3   7 


~  ^ia,p~  I  iP^ja  "  ''^ap^i7 


f     •.-     r-. 


iZ 
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Again  these  symbols  are  related  by  raising  and  lowering 
indices  in  the  usual  way  although  they  are  not  ordinary 
covariant  derivatives. 

4.  Derivation  of  the  Shell  Theory 

Let  hnlQ^jSg)  and  h2(0,,92)  t)e  negative  and  positive 
functions  respectively.   For  the  undeformed  shell  we  take  the 
region  h.   ^  0^  ^   hp-   Thus  the  shell  can  have  variable  thickness, 
and,  although  we  shall  continue  to  call  the  surface  with  position 
vector  "r  the  undeformed  middle  surface,  we  see  that  it  need  not 
be  a  middle  surface  in  the  usual  sense. 

Let  r  =  r(9,,9p)  be  the  position  vector  to  the  surface 
into  which  the  undeformed  middle  surface  goes  under  the  deforma- 
tion.  This  surface  is  called  the  deformed  middle  surface,  and 
it  is  parameterized  by  9     and  9^   in  such  a  way  that  each 
particle  has  the  same  parameter  values  on  both  surfaces. 

We  let  v"^  be  the  components  of  r  with  respect  to  the  g^'s. 
Then 

(4.1)  r  =   v^g.  =  v^g'^  . 


Letting  S^j  =  ^  a  ^®  ^^^^ 


1  1  -*/     ..I  -^x 


(4.2)  g^  =  Vl^g.  =  v.l^g 


-^        g   y  cr  .  —* 

Next  let  g^  =  -— — ^—   and  let  w"^  be  the  components  of  g^ 
with  respect  to  the  g^^'s.   Then 


^         f'.r"! 


12 


The  condition  that  g^  is  a  unit  norroal  vector  to  the 
deformed  middle  surface  yields  the  following: 

(4.4)     w^w-  =  1  ,        ^iV^I^  =  0  ,    ^-^i^^lct  =  0  . 


In  the  following  all  barred  quantities  have  the  same 
relationship  to  the  deformed  middle  surface  that  the  unbarred 
quantities  have  to  the  undeformed  middle  surface.   VJe  vjill  use 

such  quantities  without  making  further  definitions.   Indices  of 

—1 1 
such  quantities  are  raised  and  Icv/ered  using  the  quantities  g  " 

andi... 

We  consider  deformations  in  which  normals  to  the 

undeformed  middle  surface  go  into  normals  to  the  deformed  middle 

surface  with  no  change  in  length  along  normals.   This  condition 

becomes 

(4.5)  X  =  (v^+  eV)?^  . 


Next  we  calculate  P .  . .   We  have 

ox        ^   C    r  /  k  ,  n   .k  N-^  1       C(  /  k  I    ,  £,   k  I   X— »•   ,     3   k-» 


Thus 

6 


^"  hd. 


p.  .  =  (v^l  +  9  w^l  )  -^  B.}  +w^  ^  ai  . 
^ij      'a   3   'a  ox   k     dxT  k 

J  J 


Hence 


& 


.66 


15 


hi  -   -KPn.  =  <-"  I'vN„.e,wX|J  1^  a^w-^  |!2  a-} 


n         n 


^ik  J  I    'a   3   'a  cx,^      Bx   j 


n        n 


=  a"?  I  (v.  I  +e,w.  I  )   «  +  w.  -^  •  . 


n       n 


^.e, 


ot 
To  complete  the  derivation  we  compute  -r— i-  . 


X 

J 


Thus 


g^.dx  =  (6a-^3^a^^^a   ^^"^  g^ -  ax  =   de^  . 


Now  let 

(4.6)  (B^)  =  (5^-e^b^)-^  . 

Then 

d9  =  B^g^-  d5f  =  B^A^dx    and   dO  =  g^.  ^x  =  A^dx.  . 
Hence 

''^-^f  ^2  =  B^Aj;   and   ^  =  A?  . 

1  1 

Using  (4.7)  in  the  above  expression  for  P. .,  we  obtain 
(■..3)       Pia=  (^Ip+VlU'^a  .   Pl3="l- 

i      ki 

Let  Zj  =  P  P^j-  Then  from  (4.b)  and  (4.4)  we  have 


a 


z-  =    (vl^lp  +  gjwXlpjBP^   =   0    ,        z3  =   wNv^|p+ejwJp)BP  =   0   , 


•■  f 


n 


•  r 


^-I  c 


..S.1   9v    (V^V  )    L; 


e:;-i;?P 


^H 
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and 


•X,  If 

3      k 


From  this  we  see  that  the  vector  having  coir.ponents  0,  0,  1  with 
respect  to  the  g*  's  (namely  g*,)  has  a  direction  of  principal 
strain  in  the  undeformed  shell.   Since  g'^  goes  into  g,  in  the 
deformation,  it  follows  that  g^  has  a  direction  of  principal 
strain  in  the  deformed  shell,  and  the  rotation  matrix  c  rotates 
"g,  into  "g,.   In  terms  of  the  C  "^ '  s  this  means 

.'0  \    /w^'^ 

\  1  /   \  w5  / 
or 

Since  C^^C"^^  =  g^^^  =  0,  we  also  have  vj^C^^  =  0.  Kence 

(4.9)  w.C^  =  0   and   C^^  =  vj^  . 
Using  (4.5-9)  and  (3.8)  it  follows  that 

(4.10)  w.Q""""  =  0  and  Q  -^  is  proportional  to  w"^  . 

Next  we  start  the  computation  of  the  total  strain  energy 
E  =  ///  VJ(P  .,0)dv  under  arbitrary  variations  of  the  deformed 
middle  surface,  i.e.  under  arbitrary  variations  of  the  v. 's. 
We  have 

^  =/ir«''^i/^  =I/7'«'°<^ile"  Vile'^a*«""i'*"  • 


r  .-'  -r ; 


£5' 


.k'  .V   sri.-f  Id 
J. 


or. 
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Since  w.w  =  1,  we  have  w.w  =  0.   Hence  from  (4.10),  Q  -^w.  =  0 
and 

For  dv  we  have 

12;) 
Since 

If  =  S,  +  ^-,rit,   -  (  5^  9;bP  %       and   ||-  =  ?,  , 
this  becomes 

dv  =  (  6j  -  e^b^  )4  M  &^  -  93b|  )^g  .  t^d9^d9^d9^ 

=  (5^-9^bJ)(6|-e^b^)e^p^l^3.-ye^,e2de^- 


/g  (1  -  23^H+  9"K)d9,d9oda^ 

»  o  ^       ^       i    ci    ^ 


v;here  H  =  i  b^  is  the  mean  curvature  and  K  =  det  (bo)  is  the 
Gaussian  curvature  of  the  undeformed  middle  surface.  Letting 
dA  =  /g   d9nd0p  be  the  element  of  area  of  the  undeformed  middle 
surface,  we  have 

dv  =    (1  -  2e,H+  0^K)dAd9^   . 
5         3  5 

Now  define 

r  ■       "^ 

j  \ 


I  \ 


t    *      c 


16 
Then 


E  = 


(S^%J^^M^%JJdA 


SO  that  E  is  reduced  to  an  integral  over  the  undeformed  middle 
surface. 


In  the  following  we  use  the  fact  that  /g   =  /g  p  ^  . 


Then 


By  the  divergence  theorem 


s 


where  the  integration  aroung  the  boundary  is  in  the  positive 
direction  when  the  undeformed  middle  surface  is  oriented  so  that 
the  positive  side  is  the  one  to  which  g-,  points,  ds  being  arc 
length  along  the  boundary  of  the  undeformed  middle  surface. 
Let 

d  Q 
^^^•12)     ^  -  di"  '   "a  =  ^ap^   '    A  =  e^  np  . 

Then  A  g   is  a  unit  tangent  vector  to  the  undeformed  boundary 
curve,  and  n  g*   is  a  unit  outer  vector  which  is  normal  to  the 
undeformed  boundary  curve  and  is  in  the  tangent  plane  of  the 
middle  undeformed  surface. 


'US 


W*J 


•'      iJ,  ■  i. 


,{^-,  £;>;V;.    6lbf. 


;J  .  t 


;15{ 
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Thus 


1 
and 


(/gsi^J.adSi^a'/siXV^ 


Repeating  this  procedure  for  //  M  w. |  dA  we  obtain 


Next  we  observe  that  w.  is  a  function  of  the  v.]  's,  and 

1  J  "^ 


we  let 


5iW.  .  .       C;W.  .  ,       .  , 

i       ija      1       ija    ik  ja 


^""'^^^      ''i'''  ^  aTX  '   ""'"  ^  sTX  '   '^'" .  ^  ^'"''^ 


j'a  J  'a 


Then  w.  =  w.*^  v.  I   and 

1      1     J  Oi 


Methods  used  above  show  that 


M^^P|pw/-;j^d..  =7MJP|pW>X^*= 


J 

Hence 


f/(M^^l,w/-)|,V^. 


(i^.l4)   E  =/[(S^'^-M^^|pW.i°')n^:.+MJPnpW.i%J^]ds 

(si'^.MJ'^lpw/«)|^v.dA 


v-jhich  is   the  desired  form  for  E. 


■. '  r 


lisi   aw 


.    Ah,v 


sf-' 
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Let  T  g.  be  the  force  vector  applied  to  the  edge  of  the 
shell.   Then  the  principal  of  virtual  work  states  that 

where  dS  is  the  element  of  area  on  the  deformed  edge. 

Letting  s"  be  the  arc  length  along  the  boundary  curve  of 
the  deformed  middle  surface,  we  can  write  dS  in  the  form 
dS  =  Ad¥d9,  where  A  depends  on  the  v  's  and  their  first  two 


derivatives. 

Me  define 

f 

^2    _ 

(4.15) 

\ 

1 

4 

\ 

hg 

a^   =   f    9  T-^AdS 

\ 

\ 

Then  F-'""g,  and  G''''g.  are  called  the  force  vector  per  unit 
deformed  arc  length  and  the  surface  moment  vector  per  unit 
deformed  arc  length.  The  phrase  "surface  moment  vector"  is 
used  since  G"""^.  is  a  moment  aroung  the  middle  deformed  surface. 
The  phrase  "moment  vector"  is  reserved  for  the  moment  of  the 
force  vector  around  the  origin.   Then  the  principal  of  virtual 
work  becomes 

E  =  J  (F^v^+ GJwj^%.|^)d?  . 

Comparing  this  with  (4.14)  we  obtain  the  differential  equations 
(4.16)  (S^"-MJPlpW,^°^)|^  =  0 


f^rib':    r' ■''.'.■ 


^r..;. 


•ij-r  P'^ 


~      1 


,-n 


^;A^»:       ~>^ 


QiU 


«— "-^ 


*©'0i.  .. 


Siiifi- 


:».■       'IX     >• 


t! 
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and  the  boundary  conditions 


(4.17)    / 

J    ^  u 


pif  .(si-.M^P|,w/-)nJ;, 


.(G^f-M^%}w/%jJds  =  0 


for  all  variations  v. . 

iQ.      iot  i 

Since  S   and  M   are  functions  of  the  v  's  and  their 

first  two  derivatives,  we  see  that  in  (4.l6)  we  have  three 

fourth  order  equations  or  a  twelfth  order  system.   VJe  will  show 

later  that  this  system  can  be  expressed  as  a  tenth  order  system. 

Meanwhile  we  extract  what  information  we  can  from  (4.17). 

Equation  (4.17)  has  the  form 

for  arbitrary  functions  v. .   Using  the  relation 

,  .    ^  ..  .  .e 


we  have 


0  = 


[AH.  +  Bi'^(;jpA\4-v.|pn\)]ds 


5v.     .     6v, V 


V    1      a  6s        a  on  y 


where 


5v.    .     o     5v.    .     a 

5v.  _^.  .  _^. 

Then  -g—  g"""  is  the  derivative  of  v^g  "^  with  respect  to  arc 

length  in  the  direction  of  n*^^^  and  ~^   s"  is  ^  (v^'g  )  .   Hence 


:lj    b.'i  ■;, 


•-:tnr;rr;   Kr:t   af^ ; ^^      -  .'    v  '--v^: 


donoH     .  ( 
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5v. 
on  the  londeformed  boundary  curve  the  values  of  -^-^  are 

5v. 
independent  of  the  values  of  v.  and  -h— ^  .  It  follows  that 

X       OS 

B^'^n^  =  0  and 


0  = 


5v. 


A 


V,  +  B  A_  -p —  ds 
1      a  5s  y 


aH.  +4-  (B^V^J-v.  i-  (Bi°^A  )■ 
1   ds      a  1    1  6s      a 


ds 


a'  -  Is  <s'X' 


v.ds 

1 


where 


5i  (^  ^a^El  =  di  ^^  Vl^  • 


From  this  we  have  A^  =  ~  (B^^^X^).  Hence  (4.1?)  becomes 


(i^.l9) 


where 


pig=,S-.Mje|p„/«,„^.^ 


l(GJi-MJ%)w/°n„  =  0 


(GJf-MJ%)w/% 


J 


5 


[(G^f-M^-^n, 


\   lot. 
P   J    aj^i   as 


(GJg-«J%)w/Xg;J 


Equations  (4.19)  express  the  F'^'s  and  G"*" '  s  implicitly. 

Our  next  task  is  to  obtain  an  explicit  expression  for  them.   To 

i  ict 
do  this  we  derive  an  explicit  expression  for  the  W^   's  and 

several  other  relationships. 

The  relationship 


^3  = 


6l«'S2 

— ».   -»/ 


implies 


f    A^t 


)■> 


?l 
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(4.20)  w^   =  -i-  e-^'^'v,  Lv, 


g 
Then 


rz  k'l''i-'2 


-  ki         I  k.g      I         I 


and 

SI 


-ycD 


=   5^   v^l     +5     v^l       . 


3v    j  ay      -CD       oco      '7 

i  a 


Since  g  =  e„    g  ig  o>    we  have 


Sv.l  7"i       a7      '1        al      '7'°a)2     ^yl^    cx2      'oo       oo)      '2 


I'a 


=  2e     e     g     v"^! 
and 

.ii 


S        /J^\    ^  1 


2e     e     g     V    . 


1    -aoj  — V7  —     ,  i  I 
=  —  e       £       g     V 

Jg 
Thus 


Jg 
=  -^  ^V7   g^y  l^wJ    +  -i^   (^jiAl^£l2  +^jkA2\ll) 

igg 

Jgg 


■Sn^ij/ 


"•t  \' (x^jj"     S\' 


C^'--  a..c~    .C 


:^,  I 


.*  I   \ 


.^•) 


^r'T 


..).';i5 


or,-:. 


(      I    V       '^ 
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so  that 

(4.21)  wJi^  =  £Ji\|pe^P-i<^V|pwJ    .       ' 

Next 

^a  ^p  'a  1    'P^j    ijk   'a   '3^ 

On  the  other  hand 


—   —    —   l?''^   —    -»u- 


Thus 


e^^'^v.  L v.L  =  e^ow'^  . 
I'a  J ' p    ap 


Similarly 


and 


Thus 


g^xg^  =  v^i^^.  xwJj_.  =  ^^^^^^lyt^ 


Sa^S^  =  £a3pS^=  ^gaS  ^  =  ^^a^p  =  "^a^'pS 


Using  the  first  of  these  two  results  and  (4.21)  we  obtain 

Using  the  second  of  the  above  two  results  and  (4.21)  we 
obtain 


w-  Wj  =  e'pv  i^e  '  -  g  -V  Ip 


—  -CD7  i  I  •-a^  — aP  1 1 

—ay  i I   — aP  1 i    ^ 
=  g  V  I  y  -  g  ^v  I  p  =  0  . 


^f . 


«_ 


^    ^- 


J  ^\^^-     ^"'    ' 


.   0 


nl 


■:b 


airrfT 
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VJe  list  these  as 


(4.22) 


0 


wJ^%.  =  0 


Observe  that  v"*"!  A^,  v  |  n*^,  and  v;  are  the  components 
with  respect  to  the  g. 's  of  three  orthogonal  vectors.   In 

particular  v"""  |  T^g.  =  >>°^g„  is  a  unit  tangent  vector  to  the 

deformed  boundary  curve,  v  I  n^g,  =  n^g^  is  a  unit  outer  vector 

v;hlch  is  normal  to  the  deformed  boundary  curve  and  is  in  the 

— » 
tangent  plane  to  the  deformed  middle  surface,  and  w  g^  =  g,  is 

a  unit  normal  to  the  deformed  middle  surface.   Hence  there  exist 

functions  B,  C,  and  D  defined  on  the  deformed  boundary  curve 

such  that 


ds 


ds 


(4.23A)   G.  ^-M.^n^  =  B  ^  v,  |  A' +  Cv,  I  n' +  D 


T^ 


j  ds 


ds  yr 


yy 


d's 
ds 


w 


Using  (4.22)  we  obtain 


(4.23B)     (G,^-M^^.)wJi°'  = 
J  us   J   p 


-(B  -5—  A  +  Cn  )w 
ds 


Since 


d9    d0  ^ 
T-a  _  a  _  a   ds  _  ds  -.a 

ds    '^^     ds       d"s 


we  have 


TT  ^T      tP  .   g 


n,  =  e  n>^  =  —  e^o^ 


a 


'ap 


g  ^a^ 


P  ds  _  g  ds 


n. 


ds 


ds 


a 


ci. 


* ;  .t-      w  J.      *.  •>  i  «^.' 


i5    si 


..a"'.-^ 


C    j3lfS     ,0 
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We  list  these  as 

(4.24)        A°.^A-,   H„=]|i|n„. 

ds  °  ds 

Prom  (4.23B),  (4.24),  and  (4.19)  we  have 

■4  S 

Therefore  C  =  0. 
Next 

Thus  vje  can  rev</rite  the  first  of  equations   (4.19)   and 
equation   (4.23A)    as 


(4.25)       ( 


pi   =    (si^-.MjPLw>°^)n,  ^  -1-   (BwM 

''^   "^  "  ds        5s 


G^   =  n^%^^-  +  BviL A^+Dw^    . 
ds 


When  the  edge  forces  are  prescribed,  equations  (4.25) 
provide  only  four  constraints  on  the  v"""  •  s  since  the  functions  B 
and  D  are  not  determined  by  the  principle  of  virtual  work. 
Since  in  (4.l6)  we  have  a  twelfth  order  system  of  equations,  we 
are  not  yet  in  a  position  to  say  that  the  deformed  surface  is 
determined  when  the  edge  forces  are  prescribed.   The  first  step 
in  resolving  this  problem  is  to  svjitch  from  differential 
equations  and  force-deformation  relationships  expressed  with 
respect  to  the  vectors  g.  to  those  expressed  in  terms  of  the 
vectors  g".  . 


25 


5.  The  Theory  In  Terms  of  the  Deformed  Middle  Surface 

Let  F'''  and  G"""  be  the  components  of  the  force  and  surface 
moment  vectors  per  unit  deformed  arc  length  with  respect  to  the 
vectors  g.  .  Then  for  example  F'^g.  =  F"*"  g .  so  that 

F^vM^+fV  =  P^  . 

Multiplying  by  v . ] g  we  obtain 

Multiplying  by  w.  we  obtain 

F^  =  F^w.  . 


Thus 


(5.1) 


F  =  F^v.  I   ,  '^ 
a  i'a  '    a 


I  F^  =  F-^w,  ,     G3  = 


i'a 


G^w, 


Using   (5.1).    (^.24),   and   (4.25)   we  have 


a 


=    (I  v.l    (si''-MJPLw.i^)H    -v.,L  ~  (Bw^) 


I'a' 


rj     "V      i'a 


5s 


(5.2)       / 


?     =   l^w.(si^-MJ^Lw/°')K^.w,   -^   (Bw^) 


G      =    CS  V,  I    M^^p  +  B\ 


a       I—   '  1'  a         P         a 
^  g 


Equations  (5.2)  can  be  simplified  considerably  as  follows. 
From  (4.10)   w.Q^"  =  0.   Using  this  with  (4.11)  we  have 


Slic^  .o-^  ^'>?Qa.9'3 


30    ''     ~ 


jt^t;  o3v 


,8Wo/Io1 


1       1 


Prom  (^.21)  and  (4.22) 
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and 


^ay^''  ^^''  -  ^''  Sa7^'  =   ^a^'  "  ^^^'  =  °  * 


Also 


5s  ds    -^        ds  -^     >'H  X 


and 


w  4-  (Bw^)  =1-4-  (BIJ  =^  . 
5s         -^  ds    -^    ds 


With  these  results  (5.2)  becomes 


(50) 


(  F 

a 


i  g 


y-^^V^^ 


F^  = 


Ga  = 


J  g  -^  ^    ^    d¥ 


g 


^  Gj.D 


We  now  define 


J 


'  tZ''^ 


ai: 


f^i^ 


bnR 


;^^iJ5Jo:)!e?•i 


«i(     if?  r      ;  ;»        >' 


r    Si'     <i,-,; 


» 


;a  +  c 


snilein 
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'      T  f^  =  i'l  V.  I  S^^^ 
a     -  I'ct 


(5.^) 


/ 


3a 


a 


P 


The  Indices  on  the  x's  and  y-V '  s  are  to  be  raised  and 
lowered  using  the  quantities  g"'"'^  and  g^^.   Then  (5' 3)  becomes 


r 


(5.5) 


F  =   T  ^o  +  Bbg  A^ 


-3     30-   dB 
ds 

G^  =  ^W°^%  +  b"°^ 
G^  =  D 


Equations  (5.5)  are  the  desired  force-deformation  relations 
in  terms  of  the  vectors  g. .   VJe  now  start  the  calculation  of  the 
differential  equations  in  the  desired  form. 

Observe  that 


,aP 


Here  v 


a^ 


=  gr 


v'ia,.%.^M.r,^g,  =  (v^u,^rp^H,)g, 

(v-^L  )|p  as  defined  in  (3.I3).   Also 


,ap   ^a,.? 


'a"P 


'  ap  ^, 


=  rj  g^+  rj  g,  =  ( r^pv^  L + vw^  )3'^ 


Therefore 


■jHB 


9l9WOl 


(?J^) 


i^<>:'    -. 


....     .)  - 


OR  I A      .  ( .      . 


V  "■- 


'■   ^'^  ■'"1  , 
.0    <^ 


28 


(5.6) 


v'ia3  =  (r,^r>^i,.V' 


i  CL 

From  (5.4)  and  w.M   =  0  we  have 


(5.7) 


M 


^^^   '  -  v^L-'>^-^°' 


g  ^  'P 


lai 


We  define  T^^Hg  as  in  (3.13)  vjith  D^^ ,    |,  and  P^^ 
replaced  by  t"^*^,  ||  ,  and  p.^  respectively.   Since  •v'^-^  is  not 
defined,  we  define  f//   ||  to  be  an  ordinary  covariant  deriva- 
tive with  respect  to  the  metric  of  the  deformed  middle  surface, 
i.e. 


0,,  ap  I 


!'>?«? 


a  y. .  cdP  ,  ( — »  3  ^-o^^owi 


From  (5.7)  and  (5.6)  we  have  (recall  /g  ^  =  /g  n^R) 


,a 


M 


la 


S    P'^  y',a   '  JO^J  e 


r= 


p '"     "^ '  7a  jg  "  'p 


-  Jilv^lp(M^a^^,^),7^^)^  ■•>^p"(v^ip,,-r>Jip) 


f  \^H^m  ^\-  r^l  ^?r)  *  •^^^"(vMp^  .rp^vM,)] 


"  Jg  I   'P  ^   "a^°pa  ^^  w  I  . 


Vie   list  this  as 


-.sv:.-i  ii.' 
tO0S'i".;j3 


3L:  aw 


t^: 


i,.;j 
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(5.9)    m1"L  -  k  (v^U  -Wl^^iL  +  b„„  #(PV)  . 


a   J  g  ^   '  P   ^   "a   Pa 


From  (5.9)  and  (5-'^)  we  have  t   =  gc^^'^^{   L  or 
T^  =^>7^°'P||    and  from  (4.22)  and  (5.9)  we  have 


M^^l.w.i- 


-Jf'/r^llpw^ 


g 


We  list  these  as 


(5.10) 


p  J      Jg 


From  (5.4)  and  w^S'^^'  =  0  we  obtain  S^'^ 


|vM/«. 


VjfvM^ 


Hence  from  (5.10) 


Letting 


si°-MJP|pW.i"=J|(vl|pTP%T5%i)    . 


'P      j  J     g  P 


the  differential  equations  (4.l6)  become  D-^°'|   =  0,   This  means 


that 


(yg  D  Si)^a  =>^S  (TapD  g^  +  D   ^^g^  +  D   f  ^^g ^  ) 


aa 


=  yi  D-|„?,  =  0  . 


But 


\Or^     _  /— "  /,1|  .,3ct   3a  i 


/^D^^^  =ii  (v^|pT^"  +  T^"w-)f, 


g  (t   gp  +T   g^  )  =  Jg  T   g^ 


Therefore 


..,>,^ 


[SitS-iOV 


V  t 


yuK! 
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„   ,  [^  la  ^  V      f^  ia.i  — 


-^ 


let 
as  was  the  case  for  D   and  the  differential  equations 

-> 
associated  with  the  vectors  g.  are 

(5.11)  T^^^ll^  =  0  . 

icx  iot  i 

Since  S       and  M       depend  on  v     and  its  first  two  deriva- 

tives,  we  see  from  (5.4)  that  t   depends  on  v  and  Its  first 

two  derivatives  whereas  x'^     depends  on  v"""  and  its  first  three 

—  .   -*.  — *  — 

derivatives.   Since  P.^  =  g'^'g^.  ^^  vje  see  that  P.^  depends  on 

v"*"  and  its  first  tvjo  derivatives.   Therefore  (5.11)  is  a  system 

of  two  third  order  equations  and  one  fourth  order  equation 

making  (5.11)  a  tenth  order  system  of  equations.   VJe  expect  then 

that  we  need  five  boundary  conditions  to  determine  the  solution. 

After  prescribing  the  edge  forces  (5-5)  gives  us  four 

constraints.   Prescribing  the  functions  B  or  D  would  give  one 

additional  constraint  determining  the  solution  unless  buckling 

occurs.  We  next  look  for  a  reasonable  criterion  for  the  choice 

of  B  or  D.   We  v;lll  check  first  to  see  if  the  equilibrium  laws 

for  force  and  the  moment  of  force  give  a  criterion  for  such  a 

choice. 

6.  The  Equilibrium  Laws 

-i- 
Let  T  g.  be  the  edge  force  vector  per  unit  deformed  area. 

Then  the  total  edge  force  vector  is 


ij": 


ir.ei: 


^r.. 


:j   ciil^i 


(iX.?.) 


.••>ST.'f   ■?  ;";ri'j":-ni!*  ^'•.fijj    tsq 


sj;  io?09v  eoTto":  •    erl*  rscrfr 


c.ri 
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We  have 


j  T  n^g.ds  =  j  Jg  T  g^e^pAPas 


^JS   -^^""llaVV^S  =  °  • 


Also 

/(B^xPi:-ps,M3=/,-Br,?rPi-pi3M3 

as        ^      ■'^'        ds 

=  -  /  4:  (Bg^)ds  =  0  .. 
^  ds    -^ 

Thus  the  total  force  vector  is  zero  regardless  of  the 
choice  of  B. 

The  total  moment  of  force  vector  is 


T^i^x  (F  +  0^g^)dS  =J    (F^i^x  r  +  G'-i^)^  g^)ds 

=  J(Ff.x?+G%^i^)d¥ 

ds  -^        ^ 
First 


X  r 


■'-r^^ 


t''r))^   = 


tb[ 
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/.     — *         — >  r>    r—        .     — *         — ) 

T'^^'g.  X  r  H  d"s  =   /    f  g  T'^'^'g.  X  r  e     cl0 
^1'^         7  J    ^  1  ycD     CD 


=  ^(/'(  Js  T^^||^i^><  r  +  Jg  T^^i.  xiy)de3_d9. 


Next 


where  "Sap^lL  is  an  ordinary  covarlant  derivative  and  "sg-jlL  =  0* 
The  above  becomes 


u 


.0 
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/ 


Finally 


=  j    (-Bl    33^  g^^  r  -  -^  g^A  r  +  Bnpg^)ds 
=   ^    (-Ba^^  ^^^"-"=63'^^  +  Bnpg    )ds 

Q.  S 

^  ds  ^  ^ 

=  j^  (Bi^K^  +  BHp^)ds   =  f  (Bi^xi^A°'  +  Bnpi^)ds 

=  /B[e^j3A^#  +  np#3ds  =  J^El-n^g^  +  n^g^jds  =  0   . 


The  total  moment  of  force  vector  is  thus 
j      (T^^  +  b  ^^^/^^I'e^og  dA.   Since  this  is  independent  of  B  and  D, 
vje  see  that  this  will  not  yield  a  criterion  for  choosing  B  or  D. 
We  continue,  however,  since  v;e  are  interested  in  the  equilibrium 
laws  for  their  own  sake.   We  see  that  the  total  moment  of  force 
vector  is  zero  if  (t   +b  ■^''I'l      )    is  symmetric.   We  proceed  to 


show  that  this  is  the  case. 

Since  c*p  =  Jp*p,  then  c  =  piJP^)''^.      Hence  C^'^  =  P^j^D^J 
where  (D  '^)  is  symmetric.   Then  from  (3«3)  we  see  that 
Qij  ^  pi  ^  kj  ^j^gj^g  {  2^h    is  symmetric.   Since  w^Q^*^  =  0, 
w.P^^  =  0,  and  w.P^,  =  1,  then  0  =  w.Q^^'  =  w.P^,  0^^  =    ^^°'. 
Hence  Q   =  P  o  '^   where  (x   )  is  symmetric. 

Since  i^^^  =  w^l^f.  and  g^^^^  =  f  ^^I.^   =  -b^gp  =  "V^'pSi. 
then  w  I   =  -b^v  |q.   Hence  from  (4.15)  we  can  write 


P^  =  vi|p(6^.93b^)Bj.   Thus 


'.'.Uf.':'' 


f^og      e;  ' 


-.r 


;H 


3* 


r"P=r''T/  =  g="''  f  vj/P 


V  g 


Let 


.j-   7  -    o)        3     5 


Then  (D   )  Is  symmetric  and 


Similarly 

Hence 

7'     J  57      37      37Cd': 

=  r(D«P-92^XV)d9,  . 

Since  (D   )  and  (bTD^'^V^^)  are  symmetric,  the  desired 
conclusion  follows.   Hence  the  total  moment  of  force  vector 
is  zero  and  the  usual  equilibrium  lav;s  are  satisfied. 


f.» 


^   ,.1  -: 
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7.  Choice  of  B 

VJe  finally  obtain  a  reasonable  criterion  for  choosing  B  by 
taking  a  close  look  at  the  parts  of  the  force  vector  and  surface 
moment  vector  per  unit  deformed  arc  length  which  do  not  depend 
on  B  or  D.   First  from  (5.'0 

P-    fg    I  ciP-      I  qIP   ds 

— » 

Let  N  be  the  unit  outer  normal  vector  and  dS  the  element 
of  area  for  the  undeformed  edge.   Then 

ifdS  =  f-  (r^e^gy  x^  {r\  9^g^.)dsde^ 


Let  N°^„  =  t.      Then 


=a 


N^dS  =    (5^.93b^)A^'e^^dsd93 


Multiplying  by  5'  -  9^b^  we  obtain 

jr      a  a  u  ^  ^     (Xi 

=  /ie^det    (5j-e3bJ)dsd93 

=   n   (1  -  2H9-,  +  K95)dsd9,    . 
CO  3  3  3 


Writing  dS  as  AdsdQ.,  v;e  have 

3 


!    ..T 


OViJri    9V/ 


(7.1)  N^(5^-S^i^^   =   nJl-2H93  +  K9|)    . 

We  now  have 

dS 
Recall  from  (2.2)  that  q. .n.  —  are  the  Cartesian 

^^  J  dS 
coordinates  for  the  edge  force  given  by  the  three  dimensional 

theory.   In  Section  2,  (n  ,np,n^)  was  the  unit  outer  normal 

vector  so  the  n . ' s  there  are  different  from  the  n  's.   Hence 

p,ia„     dS  -»  aIo'tvt     Ads  -*  ,  ,-^       ir^—^a    \ 

^     N     -—  g     =   Q     N     ::: —  g  (dS  =  Adsd9,) 

^  dS     ^  ^  Ads     ^  ^ 

is  the  edge  of  force  vector  as  computed  by  the  three  dimensional 

theory.  Let  Tn   =  ^n  U^^^No  ^  and  "T^   =  w,  Q^^No  ^.  Then 
-  ^^  ^^    ^  "^    ^  Ads        ^    ^    P  Ads 

^     g.  is  the  edge  force  vector  as  computed  from  the  three 

dimensional  theory.   Since  w.Q''"°'  =  0,  then  'S       -   0   and 

—2  i  —    ,-~-<x  — 

\r  Sj_  =  J   S^j  is  tangent  to  the  deformed  middle  surface. 

Then  t  ^n^  =  /  TId9,.   Hence  r^'^no  =  /  T°^de_  are  the 
apj-'a5  Pj       3 

components  of  force  tangent  to  the  deformed  middle  surface  per 
unit  deformed  arc  length  as  computed  from  the  three  dimensional 
theory.   Hence  if  we  choose  B  =  0  in  (5-5),  the  stress  strain 
laws  of  the  three  dimensional  theory?-  hold  in  the  direction 
tangent  to  the  deformed  middle  surface.   VJe  make  that  choice 
so  that  {5. 5)  becomes 


(I.V) 


ta   9f. 


rto2;jo9'ilb 


-  r 
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C  m    ia— 


(7.2)  <  G^  =  %f% 


^ 


i^  G^  =  D  . 


Then  it  also  folloivs  as  above  that  the  surface  moment 

vector  as  computed  from  the  three  dimensional  theory  is  tangent 

— » 
to  the  deformed  middle  surface  and  is  G  g  .   The  component  of 

force  normal  to  the  deformed  middle  surface  is  not  consistent 

with  the  stress  strain  laws  of  the  three  dimensional  theory, 

but  it  is  such  that  the  equilibrium  laws  are  satisfied.   The 

component  of  the  surface  moment  vector  normal  to  the  deformed 

middle  surface  is  also  in  general  not  consistent  with  the  stress 

strain  laws  of  the  three  dimensional  theory.   Its  value  is  given 

algebraically  in  terms  of  the  prescribed  edge  force  and  the 

deformation  by  G-^  =  G"^w.  . 

8.  Comparison  with  the  von  Karman-F8ppl  Equations 

Consider  a  homogeneous  isotropic  plate  vjhich  occupies  the 

region  |x,i  _^  -p  in  the  undeformed  state  so  that  h  is  the 

thickness.   Call  x-,  =  0  the  middle  undeformed  surface. 

3 

Choose   9     =  X   .      Then  6     =  x,  from   (3.2)   and  k^.   =  a'^   =  b\ 
"'  o,  y  t>  J  J  J 

Hence   P..   =  p.-,    C.   =  c..,    Q.  .   =  q.-,    and   indices,    vjhich  are 

raised  and  lov;ered  v/ith  the  g"""'-    and  g.  .,    have   the   same  meaning 
v;hen  they  are   up  as   v;hen  they  are  down.      We  will   write    such 
indices   in  the   lower  position. 


ns" 


Siii."rycn'.   ^xpi.   5d:J   svsri 


5u 


Since  g^^  =  ^^ii'^i2'^i3^^  ^^°"^  (^.1),  (^-3),  and  (3.10) 

-  -  t-i  J 

we  have  r  =  ( v-j^,  v^,  v^) ,  g,  =  (w^jW^^w^),  and  \    .^^  =   0.   Hence 

V. I   =  V.   ,  etc. 
I'a    i,ct 

The  following  are  obtained  from  (4.11),  {K.lC),    (4.20), 

and  (4.21): 

h/2 


(8.1) 


/ 


la 


-h/2 


q.  dx-,  ,   M.   =  /    x,q.  dx. 


h/2 

r 

-h/2 


(S.   -M.o  oW..  )    =  0 


■>  Sap  =  ^i,a^i,3 


w. 


e.  .,  V  .  T  V,  ^ 
ijk  j,l  k,2 


^  ^jik\,P^ap-^S,P^j 


In  the  usual  derivation  of  the  von  Karman-F'Bppl  equations, 

in  addition  to  the  Kirchhoff  hypothesis  made  already,  it  is 

2 
assumed  that  v  -  x  and  their  derivatives  are  of  order  h  and 

ct  a. 

V-,  and  its  derivatives  are  of  order  h.   VJe  incorporate  these 
3 

assumptions  by  assuming  v.  is  asymptotic  to  a  power  series  in  h. 


00 


^=1 


n=0 


V.  h 
m 


n 


v/here  v   =x,0=v,=v^  and  that  differentiating  these 
ao    a'  al    ^o 

series  term  by  term  with  respect  to  x^  yields  a  series  v/hich 
is  asymptotic  to  the  derivatives  of  v,.  . 


For  convenience  we  let  v  =  v 


31' 


Then  from  (8.1) 


r]&i.'.\v 


n 


a> 
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{^ 


Sap  =   ^ap+(^a2,p+ V,a''\a\p)^    +  O^^''^ 


P 


-1 


(S^     )     =     (g„;^)"        =     6„R-    (V^o    R+V«o    „+V    ^V    ^)h      +0(h^) 


=ai3 


g  =   1  +  (2v^„   ^  +v     V   ^)h    +  O(h^) 


ap      ^    a2,3"^32,a"\a\P' 


1^2  J-  n^v^5' 


J. 


72,7         ,7    ,7 


'-^V,74\7\7^^    -^^^^^ 


g 


(8.2) 


^a  = 


w^  = 


•\a'^-  V,a^^+0(h^') 


1-^  h^v      V   ^  +  0(h^) 


^Pra  =  °^^^ 


W3pa  =   0(h    ) 


^P3a=   -^aP+O^h-) 


i.  ^33a  =  °(^) 


From   (4.0)      p.^  =  ^i,a  + Vi,a  ^^^^  ^13  =  ^    ^°   ^hat 


^ag 


^ap-^3^\ap+h\2,p  +  °(h^) 


I   P3a 


(6.3)     <j 


^a3 


V  ^h+ V  h"  +0(h5) 


-\ah-^32,a^'-^0(h5) 


,  P33=  l-|h2v^^v^^  +  0(h5) 


^^^  "ij   =  PkAj  -  ^ir      Then  z.^  =  z^.    =  0,   and 
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(8.4)   z^p  =  0(h2)  =  -2x3hv^^p  +  h2(v^^2,^-^^^2,a  +  ^a^^)+0(h5) 

Let  p,  Cy    2,  and  q  be  the  matrices  (p.  J,  (c..),  (z,.), 

-'-J     -*-  J     J- J 

and  (q,-,-)  respectively.   Then 


(8.5) 


-B 


c"p  =  ip*p  =  /r+T  =  1+1  z  +  o(h""'") 


s-^  =  [c*p-l]  =  I  [z]  +0(h^) 


,2i    ^,JU  _    ^/,6. 


32  =  [{c*p  -1)^]  =  0(hM  ,    s^  =  0(h°) 


c  =  p(j'p^)"^  =  p[l-|  z  +  0(h^)] 


where  the  square  bracket,  [   ],  denotes  the  trace  of  the  matrix. 

For  a  homogeneous  isotropic  material  we  want  the  strain 
energy  density  W  to  agree  to  lowest  order  terms  with  the  strain 
energy  density  of  the  linear  theory.   That  is,  we  want  to  have 


A   2 


V/  =  -^  s,  +11S2  +  higher  order  terms  in  s,  ,  Sp,  s. 


where  >v  and  n  are  the  Lame  constants.   It  follows  from  (2.7) 
that 

q  =  W^p  +  UgC  +  V;  pc*p 

=  p  [ll^  +\^^[1  -\   z  +  0(h'^)]  +W^[1  +  -|  z  +  0(h^)]  J 

=  p  [-hs^  +nz  +  0(h'h]]. 


This  with  (G.l)  gives 


••^) 


.^) 


\  \C\l)C  h- 
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Z' 


(8.6) 


^a3 


^3a 


'ap 


=   0(11^) 


h^[A5^.(v^,o  .  +4  V     V   ^ J 


^ 


M 


aP 


3a 


M 


aP^>2,7     2    \7    ,7' 
+  ^(^a2,P  +  ^P2,a+\a\p)^+0(hS 

-fe   (^V^^+2^tv^^p)  +  0(h5) 

O(h^) 


=   O(h^) 


Where  Av  =  v 


J  77 
Defining  A  g  as  below,    (8.6)    and    (o.2)   yield 


(8.7) 


A 


ap 


s     ai3 


iV 


'^5ap(V,7^2   \7\7^-'^^'V,P+^f32,a-^\a\3) 
S^p+0(h5)    =   h\p+0(h^) 


Hence  A  g  r  ~  '-'  ^'^'^  since  A  g  =  Ap  ,  there  is  a  function 


(|)   such  that 


^11  ^ 


(8.8)^    i 


,22 


^22   =   ^(^12,l-*-|  \l\l^+  (A+2n)(v22^2  +  l  ^  2^  2  ^ 
^11      =    (^+2;i)(v^2, 1  +  1  \l\l)+  ^(^22,2+1  ^2^2^ 


-^12  ^   ^12  =  '^^^12,2+^22,l+^1^2) 


From  the  first  two  of  (8.8) 


n 


^\"'- 


^a.Q^^^-*-x)f- 


's.^/-*-;,. 


k2 


-l2,l4^l^l  =  " — — — ^ 


^22,2-'2  ^,2^,2 


(A+2:.l)(1)  ^^  -  A(l)  22 
~    ^l-H  (  A+y. ) 


Prom  these  and  the  last  of  (3.8) 


■''',1122   ^12,12 


=  f^(^2, 122  •'^22, 112 +^,112^,2  +  ^,11^,22 


+^l2^l2+^l^l22) 


=  \^ 


'(^+^^^^ 'I',  2222-^^1122  .  ^  ^+^^^  ^ ^,  nil -^'>\  1122 


-    \l2\l2+\ll\22 


(8.9)  A2^=^^iiliil(v^^2\l2-\ll\22^ 


E 


1  -cr 


2    ^\12^,12-^,11^,22) 


v/here  E  is  Young's  modulus,    J' is  Poisson's   ratio  and 


A^cl)   =   (f) 


,aapp* 


From  before    (U.6)   q  =  p(  As,  +  |az )  +  0(h    )    so   that 


^3a  =  hv^Q(As^5g^+nz^^)  +  0(h    ) 


=   h\p  'lh\J^-^5'^('^V-^^+^^^a^)^  +  °^h')    > 


and 


Sja  =  h\3^p+0(h5)    . 


J  add   C'S   {   ii)0  •  ' 


t    {%'.)0  ^l  \ 


.    { 


From  this  and    (G.7),    (0.6) 


^5 


N^a  =  h\pV^p  +  M^p^p  +  0(h5) 


^'^V^p-^(^^),al-0(h5)    , 


and 


yields 


0  =   N 


3a,a 


h+^  a2„  _ 


^^.,\,Ka  -  ^W  ^^^  -  ^   ' 


(8.10)      A' 


V    = 


12 


A+2iI   ('^11^22  +^22^11  -  2^1^,12^,12) 
12(l  +  '7')(i_  2r) 


E(1-t) 


(^11^22+^22^,ll-2'l^,12^,12) 


Let 


T, 


ap  -  I TJT  "^aiPgij 


X3=0 


Then  (see  (2,5))  the  T  g  are  the  components  of  the  stress  tensor 
on  the  middle  surface  as  computed  from  the  three  dimensional 
theory.   We  obtain 


'a|3  "  ^ap 


x,=0 


+  0(h^)  -  h2A^p+  O(h^)  . 


Let  V  =  h  (j).   Then  ^  is  the  usual  stress  functi 
satisfying  T^^  =  V'  22+ ^^^^^^  ^22  "  ^  ]_3^+0(h^), 

m     _  m 
^12     ^21 


on 


-f   ^2  i-O(h^) 


w 


Let  w  =  V  .  Then  v  =  ^+0(h)  and  (8.9)  and  (o.lO)  become 


onB 


•f 


i'id 
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(  .2 


(8.11)  1 


Equations  (8.11)  differ  from  the  von  Karman-F8ppl  equations 

in  the  constants  tt   and  — ^ — '—^ — ~ — —   which  are  E  and 

1  -il  h^E(l  -  r) 

^^"  ' — -  respectively  for  the   von  Karrnan-Pttppl   equations. 
h'^E 


i:r:B  : 
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ENG-EB  Port.  Constr.  Br., 

Eng. Div. Mil. Cons,  (l) 
ENG-EA  Struc.  Br.  Eng. 

Div.  Mil.  Constr.  (1) 
ENG-NB  Spec.  Engr.  Br., 

Eng.  R/D  Div.      (1) 

CO.  Engin.  Res.  Dev.  Lab. 
Fort  Bel voir,  Virginia 


(1) 


Office  of  Chief  of  Ordnance 

Dept.  of  Army,  Washington  25,  D.C. 

Attn:   Res.  and  Mat.  Br., 

(Ord.  R/D  Div. )  (1) 

Office  of  Chief  Signal  Officer 
Dept.  of  Army,  V.'ashington  25,  D.C. 
Attn:   Engin.  and  Techn.  Div.    (1) 

CO,  VJatertown  Arsenal 
Watertown,  Massachusetts 
Attn:   Lab.  Div. 


CO,  Frankford  Arsenal 
Bridesburg  Station 
Philadelphia  37,  Penna. 
Attn:   Lab.  Div. 

Office  of  Ordnance  Research 
2127  Myrtle  Dr.,  Duke  Station 
Durham,  North  Carolina 
Attn:   Div.  of  Engin.  Sci. 

CO,  Squier  Signal  Lab. 
Fort  Monmouth,  N.  J. 
Attn:   Comp.  and  Mat.  Br. 


(1) 


(1) 


(1) 


(1) 


Chief  of  Naval  Operations 

Dept.  of  Navy,  Washington  25,  D.C. 

Attn:   Op  91  (1) 

Op  03EG  (1) 

Commandant,  Marine  Corps 

Headquarters,  USMC 

Washington  25,  D.  C.  (1) 
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Chief, 

Bureau  of  Ships 

Dept. 

of  Navy 

Washington  25,  D.  C. 

Attn: 

Code  106 

(1) 

Code  312 

(5) 

Code  320 

(1) 

Code  370 

(1) 

Code  375 

(1) 

Code  420 

(1) 

Code  421 

(1) 

Code  423 

(2) 

Code  425 

(1) 

Code  440 

(1) 

Code  442 

(2) 

Code  443 

(1) 

Code  ^25 

(1) 

Code  633 

(1) 

Chief, 

,  Bureau  of  Aeronautics 

Dept. 

of  Na-.-y 

Washington  25,  D.  C. 

Attn: 

AE-4 

(1) 

AV-34 

(1) 

AD 

(1) 

AD-2 

(1) 

RS-7 

(1) 

RS-8 

(1) 

SI 

(1) 

AER-126 

(1) 

Chief, 

,  Bureau  of  Ordnance 

Dept. 

of  Navy 

Washington  25,  D.  C. 

Attn: 

Ad3 

(1) 

Re 

(1) 

ReS 

(1) 

ReU 

(1) 

ReS5 

(1) 

ReSl 

(1) 

Ren 

(1) 

Spec. 

Proj.  Office,  Bur.  Ord 

^ 

Dept. 

of  Navy 

Washington  25,  D.  C. 

Attn: 

Missile  Br, 

(2) 

Chief, 

,  Bur.  Yards  and  Docks 

Dept. 

of  Navy 

Washii 

igton  25,  D.  C. 

Attn: 

Code  D-202 

(1) 

Code  D-202. 3 

(1) 

Code  220 

(1) 

Code  D-222 

(1) 

Code  D-410C 

(1) 

Code  D-440 

(1) 

Code  D-500 

(1) 

Commanding  Officer  and  Director 

David  Taylor  Model  Basin 

Washington  7,  D.  C. 

Attn:   Code  l40  (l) 

Code  600  (1) 

Cede  700  (1) 

Code  720  (1) 

Code  725  (l) 

Code  731  (1) 

Code  740  (2) 

CO,  U.S.  Naval  Ordnance  Lab. 
White  Oak,  Maryland 
Attn:   Techn.  Library  (2) 

Techn.  Eval.  Dept.       (1) 

Director,  Materials  Lab. 

N.Y.  Naval  Ship:yard 

Brooklyn  1,  N.  Y.  (1) 

CO,  Portsmouth  Naval  Shipyard 
Portsmouth,  New  Hampshire       (2) 

CO,  Mare  Island  Nav.  Shipyard 
Vallejo,  California  (2) 

CO  and  Director 

U.S.  Nav.  Electron.  Lab. 

San  Diego  52,  California        (1) 

Officer-in-Charge 
Nav.  Civ.  Engin.  Res. 
and  Eval.  Lab. 
U.S.  Nav.  Constr.  Battal.  Center 
Port  Hueneme,  California        (2) 

Dir.,  Nav.  Air  Experimental  Sta. 
Nav.  Air  Mat.  Center,  Nav.  Base 
Philadelphia  12,  Penna. 
Attn:   Materials  Lab.  (1) 

Structures  Lab.  (1) 

Officer-in-Charge 

Underwater  Explos.  Res.  Div. 

Norfolk  Naval  Shipyard 

Portsmouth,  Virginia 

Attn:   Dr.  A.  H.  Keil  (2) 

CO,  U.S.  Nav.  Proving  Ground 
Dahlgren,  Virginia  (1) 

Supr.  of  Shipbuilding 

USN  and  Nav.  Inspec.  of  Ordnance 

General  Dynamics  Corp. , 

Electr.  Boat  Div. 

Groton,  Connecticut  (1) 

Supr.  of  Shipbuilding 

USN  and  Nav.  Inspec.  of  Ordnance 

Newport  News  Shipbuilding 

and  Dry  Dock  Co. 
Newport  News,  Virginia         (1) 
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Supr.  of  Shipbuilding 
USN  and  Nav.  Inspec.  of  Ordnance 
Ingalls  Shipbuilding  Corp. 
Pascagoula,  Mississippi      (1) 

CO,  U.S.  Nav.  Admin.  Unit 

MIT,  Cambridge  39,  Mass.     (1) 

Officer-in-Charge 
Postgrad.  School  for 

Naval  Officers 
Webb  Inst,  of  Nav,  Arch. 
Crescent  Beach  Rd. 
Glen  Cove,  L.I. ,  N.Y.        (l) 


Supt.,  Nav.  Gun  Factory 
Washington  25,  D.  C. 


(1) 


Comm. ,  Nav.  Ordnance  Test  Sta. 
China  Lake,  California 
Attn:   Physics  Div.  (l) 

Mechanics  Div.        (1) 

CO,  Nav.  Ordnance  Test  Sta. 

Underwater  Ordnance  Div. 

3202  E.  Foothill  Blvd. 

Pasadena  8,  California 

Attn:   Struc.  Div.  (1) 

CO  and  Director 

U.S.  Nav.  Engin.  Exp.  Station 

Annapolis,  Maryland  (1) 

Supt.  U.S.  Nav.  Postgrad.  School 


Monterey,  California 
Comm.  Marine  Corps  Schools 


(1) 


Quantico,  Virginia 
Attn:   Dir. ,  MC  Dev, 


Center   (1) 


Comm.  Gen. ,  USAF 
Washington  25,  D.  C. 
Attn:   Res.  and  Dev.  Div. 

CO,  Air  Material  Command 
Wright-Patterson  APB,  Ohio 
Attn:   MCREX-B 

Structures  Div. 

CO,  USAF  Inst,  of  Technology 
Wright-Patterson  APB,  Ohio 
Attn:   Chief,  Appl.  Mech. 
Group 

Director  of  Intelligence 

Headquarters,  USAF 

Washington  25,  D.  C. 

Attn:   PV  Br.  (Air  Targ.  Div)(l) 

CO,  A.F.  Office  Sci.  Research 

Washington  25,  D.  C. 

Attn:   Mechanics  Div.        (l) 


(1) 


(1) 
(1) 


(1) 


U.S.  Atomic  Energy  Commission 
Washington  25,  D.  C. 
Attn:   Dir.  of  Research 

Dir. ,  Nat.  Bur.  of  Standards 
Washington  25,  D.  C. 
Attn:   Div.  of  Mechanics 

Engin.  Mech.  Sect. 

Aircraft  Structures 

Comm.,  U.S.  Coast  Guard 

1300  E  St. ,  NW 

Washington  25,  D.  C. 

Attn:   Chief,  Test  and  Dev.  Div 

U.S.  Maritime  Administration 
General  Admin.  Office  Bldg. 
441  G  St.,  NW 
Washington  25,  D.  C. 
Attn:   Chief,  Div.  Prelim. 
Design 

Nat.  Aero,  and  Space  Admin. 
Langley  Research  Center 
Langley  Field,  Virginia 
Attn:   Structures  Div. 

Nat.  Aero,  and  Space  Admin. 

1512  H  St. ,  NW 

Washington  25,  D.  C. 

Attn:   Loads  and  Struc.  Div. 

Director,  Forest  Prod.  Lab. 
Madison,  Wisconsin 

Federal  Aviation  Agency 
Dept.  of  Commerce 
Washington  25,  D.  C. 
Attn:   Chief,  Air  Engin.  Div. 

Chief,  Air.  and 

Equip.  Div. 

National  Science  Foundation 


(2) 


(1) 
(1) 

(1) 


(1) 


1520  H  St, 


m-J 


Washington,  D.  C. 

National  /Academy  of  Sciences 

2101  Constitution  Ave., 

Washington  25,  D.  C. 

Attn:   Dir.  ,  Com.m.  on  Ships 
Struc.  Design 
Exec.  Secy.,  Comm.  on 
Undersea  Warfare 

General  Dynamics  Corp. 
Electr.  Boat  Div. 
Groton,  Connecticut 

Newport  News  Shipbuilding  and 
Dry  Dock  Co. 

Newport  News,  Virginia 


(1) 

(2) 

(2) 
(1) 

(1) 
(1) 

(1) 
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Ingalls  Shipbuilding  Corp. 
Pascagoula,  Mississippi      (1) 

Prof.  Lynn  S.  Beedle 

Fritz  Engineering  Lab. 

Lehigh  University 

Bethlehem,  Penna.  (1) 

Prof.  R.  L.  Bisplinghoff 
Dept.  of  Aero.  Engineering 
Massachusetts  Inst,  of  Techn. 
Cambridge  39,  Massachusetts   (1) 

Prof.  H.  H.  Bleich 

Dept.  of  Civ.  Engineering 

Columbia  University 

New  York  2?,  N.  Y.  (1) 

Prof.  B.  A.  Boley 

Dept.  of  Civ.  Engineering 

Columbia  University 

New  York  27,  N.  Y.  (1) 

Dr.  John  F.  Brahtz 

Southern  California  Labs. 

Stanford  Research  Institute 

820  Mission  St. 

South  Pasadena,  California   (1) 

Dr.  D.  0.  Brush 

Struc.  Dept.  53-13 

Lockheed  Aircraft  Corp. 

Missile  Syst.  Dlv. 

Sunnyvale,  California        (l) 

Prof.  B.  Budiansky 

Dept.  of  Mech.  Engineering 

School  Appl.  Sciences 

Harvard  University 

Cambridge  38,  Massachusetts   (1) 

Prof.  Herbert  Dereslewicz 

Dept.  of  Civ.  Engineering 

Columbia  University 

632  W.  125th  St. 

New  York  27,  N.  Y.  (1) 

Prof.  D.  C.  Drucker,  Chmn. 

Dlv.  of  Engineering 

Brown  University 

Providence  12,  Rhode  Island   (1) 

Prof.  John  Duberg 

Dept.  of  Civ.  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Prof.  J.  Erlcksen 
Mech.  Engineering  Dept. 
Johns  Hopkins  University 
Baltimore  18,  Maryland       (1) 


Prof.  A.  'C.  Eringen 

Dept.  Aero.  Engineering 

Purdue  University 

Lafayette,  Indiana  (l) 

Prof.  W.  Plugge 

Dept.  of  Mech.  Engineering 

Stanford,  California  (1) 

Mr.  M.  Goland,  VP  and  Dlr. 
Southwest  Research  Institute 
8500  Culebra  Rd. 
San  Antonio  6,  Texas  (1) 

Prof.  J.  N.  Goodler 

Dept.  of  Mech.  Engineering 

Stanford  University 

Stanford,  California  (1) 

Prof.  L.  E.  Goodman 
Engineering  Experimental  Sta. 
University  of  Minnesota 
Minneapolis,   Minnesota         (l) 

Prof.  M.  Hetenyl 

The  Technical  Institute 

Northwestern  University 

Evanston,  Illinois  (1) 

Prof.  P.  G.  Hodge 

Dept.  of  Mechanics 

Illinois  Inst,  of  Technology 

Chicago  16,  Illinois  (l) 

Prof.  N.  J.  Hoff,  Head 

Dlv.  Aeronautical  Engineering 

Stanford  University 

Stanford,  California  (1) 

Prof.  W.  H.  Hoppmann,  II 

Dept.  of  Mechanics 

Rensselaer  Polytechnic  Inst. 

Troy,  New  York  (1) 

Prof.  Bruce  G.  Johnston 

University  of  Michigan 

Ann  Arbor.  Michigan  (l) 

Prof.  J.  Kempner 
Dept.  of  Aero,  Engineering 
and  Appl.  Mechanics 
Polytechnic  Inst,  of  Brooklyn 
333  Jay  St. 
Brooklyn  1,  N.  Y.  (1) 

Prof.  H.  L.  Langhaar 
Dept.  of  Theoretical  and 
Applied  Mechanics 
University  of  Illinois 
Urbana,  Illinois  (1) 
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Prof.  B.  J.  Lazan,  Dj rector 
Engineering  Experimental  Sts. 
University  of  Minnesota 
Minneapolis  1'+,  Minnesota 

Prof.  E.  H.  Lee 
Div.  of  Appl.  Mathematics 
Brovm  University 
Providence  12,  Rhode  Island 

Prof.  George  H.  Lee,  Dir.  of 
Rensselaer  Polytechnic  Inst. 
Troy,  New  York 

Mr.  S.  Levy 

GE  Electr.  Research  Lab. 
5198  Chestnut  St. 
Philadelphia,  Penna. 

Prof.  Paul  Lieber 
Geologj'  Department 
University  of  California 
Berkeley  h,    Calif oi^nia 

Prof.  Joseph  Marin,  Head 
Dept.  Engineering  Mechanics 
College  of  Engin.  and  Arch. 
Pennsylvania  State  University 


(1) 


(1) 
Res, 

(1) 


(1) 


(1) 


University  Park,  Penna. 

Prof.  R.  D.  Mindlin 
Dept.  of  Civ.  Engineering 
Columbia  University 

632  w.  125th  St. 

New  York  2?,  N.  Y. 

Prof.  Paul  M.  Naghdi 
Building  T-7 
College  of  Engineering 
University  of  California 
Berkeley  4,  California 


(1) 


(1) 


(1) 


William  A, 


Nash 


Prof 

Dept.  of  Engineering  Mechanics 
University  of  Florida 
Gainesville,  Florida         (1) 

Prof.  N.  M.  Ne\^mlark,  Head 
Dept.  of  Civ.  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Prof.  E.  Orowan 
Dept.  of  Mech.  Engineering 
Massachusetts  Institute  of  Techn. 
Cambridge  59,  Massachusetts   (1) 

Prof.  Aris  Phillips 

Dept.  of  Civ.  Engineering 

15  Prospect  St. 

Yale  University 

New  Haven,  Connecticut       (1) 


Pi-of.  VJ.  Frager,  Chmn. 

Phys.  Sci.  Council 

Brown  University 

Providence  12,  Rhode  Island     (1) 

Prof.  J.  R.  M.  Radok 
Dept.  of  Aero  Engineering 

and  Appl.  Mechanics 
Polytechnic  Inst,  of  Brooklyn 
335  Jay  St. 
Brooklyn  1,  N.  Y.  (1) 

Prof.  E.  L.  Reiss 

Inst,    of  Mathematical   Sciences 

New  York  University 

•I  VJashington  Place 

New  York  3,  N.  Y.  (1) 

Prof.  E.  Reissner 
Dept.  of  Mathematics 
Massachusetts  Inst,  of  Technology 


Cambridge  39,  Massachusetts 

Prof.  M.  A.  Sadowsky 
Dept,  of  Mechanics 
Rennselaer  Polytechnic  Inst. 
Troy,  New  York 

Prof.  B.  W.  Shaffer 
Dept.  of  Mech.  Engineering 
New  York  University 
University  Heights 


(1) 


(1) 


New  York 


N.  Y, 


Prof.  J.  Stallmeyer 
Dept.  of  Civ.  Engineering 
University  of  Illinois 
Urbana,  Illinois 

Prof.  Eli  Sternberg 
Dept.  of  Mechanics 
Brown  University 
Providence  12,  Rhode  Island 

Pi^of.  T.  Y.  Thomas 
Grad.  Inst.  Math,  and  Mech. 
Indiana  University 
Bloomington,  Indiana 

Prof.  S.  P.  Tim.oshenko 
School  of  Engineering 
Stanford  University 
Stanford,  California 

Prof.  A.  S.  Velestos 
Dept.  of  Civ.  Engineering 
University  of  Illinois 
Urbana,  Illinois 


(1) 


(1) 


(1) 


(1) 


(1) 


(1) 
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Dr.  E.  Wenk 

Southwest  Research  Institute 

8500  Culebra  Rd. 

San  Antonio,  Texas  (1) 

Prof.  Dana  Young 

Yale  University 

New  Haven,  Connecticut       (l) 

Prof.  R.  A.  Di  Taranto 

Dept.  of  Mech.  Engineering 

Drexel  Institute 

32nd  and  Chestnut  Streets 

Philadelphia,  Penna.         (1) 

Mr.  H.  K.  Koopman,  Secy. 
Welding  Res.  Council 

Engineering  Foundation 
29  w.  39th  St. 

New  York  l3,  N.  Y.  (2) 

Prof.  Walter  T.  Daniels 
School  of  Engin.  and  Archit. 
Howard  University 
Washington  1,  D.  C.  (1) 

Conun. ,  (Code  755) 

U.S.  Naval  Ordnance  Test  Sta. 

China  Lake,  California 

Attn:   Techn.  Library        (1) 


Prof.  J.  E.  Cermak 
Dept.  of  Civ.  Engineering 
Colorado  State  University 
Port  Collins,  Colorado 

Prof.  W.  J.  Hall 
Dept.  of  Civ.  Engineering 
University  of  Illinois 
Urbana,  Illinois 

Dr.  Hyman  Serbin 
Design  Integration  Dept. 
Hughes  Aircraft  Co. 
Culver  City,  California 


Comaancer,  ■.'iALli 
Wrlght-Pa-iterson  AFB, 
Attn :  '■j'wRC 


Cfhio 


(1) 


(1) 


(1) 


(1) 
(1) 
(1) 


Legislative  Reference   Service 
Library  of  Congress 
VJashington  25,   D.    C. 
Attn:      Dr.    E.    Wenk 


(1) 


Dr.  A.  Ross 

Aircraft  Nuclear  Propulsion  Dept. 

General  Electric  Co. 

Cincinnati  I5,  Ohio  (1) 

Dr.  F.  Lane 

General  Applied  Science  Labs. 
Stewart  and  Merrick  Avenues 
Westbury,  L.I.,  N.  Y.  (1) 

Commander  Edward  Leonard 

Asst.  Navy  Representative 

MIT  Lincoln  Lab. 

Lexington  73>    Massachusetts  (1) 

Prof.  '.;,   Pohle 
Grad,   riath.   Depto 

Adelphi   College 

Garden  City,   L.I,,   N.Y,  (l) 

Pr,  wartin  GRldbcrg 

Research  Dept. 

Grumann  Aircraft 

Bethpage,  Ld,. ,  K,Y,  (l) 

Prof.  S.  D.  Larmarajan 

Aero   Space  Dept, 

3an  Diego  State  University 

San  Diego  15,  Calif orriia  (l) 

Professor  R,  P.  Harrington,  Head 
Dept.  of  Aeronautical  Engineering 
University  of  Cincinnati 
Cincinnati  21,  Ohio  (1) 

Professor  Eugene  J.  Brunelle,  Jr. 
Dept,  of  Aeronautical  Engineering 
Princeton  University 
Princeton,  New  Jersey 


Commanding  Officer 

USNNOEU 

Klrtlar.d  Air  Force  Base 

Albuquerque,  New  Mexico 

Attn:   Code  20 

(Dr.  J.N.  Brennan) 


(1) 


:  \ 


^^^^ 


v^ 


\# 


DATE  DUE 

GAVLORD 

PRINTEOIN  u    S   A- 

(I) 


"■'uc 


""^S'ir- . 


JV, 


v^ 


